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We show that lattice results on four flavor QCD at nonzero temperature and baryon 
density compare well with the hadron resonance gas model up to T ~ 0.95T C , and 
approach a free field behaviour with a reduced effective number of flavor for T > 
Oh 1-5T c ; chiral symmetry and confinement are interrelated, and we note analogies 

between the critical line of QCD and that of simple models with the same global 
i-O ' symmetries. 



j_] . 1. QCD Thermodynamics and Imaginary fis 

c3 ' 

In the last four years a few lattice techniques proven successful in QCD 
thermodynamics for hb/T < 1. 1 ' 2 - 3 < 4 ' 5 ^' 7 . While waiting for final results in 
the scaling limit and with physical values of the parameters, it is very useful 
to contrast and compare current lattice results with model calculations and 
perturbative studies. 

The imaginary chemical potential approach 8 ' 9, 5 ' 1,2 ' 10 to QCD thermo- 
dynamics seems to be ideally suited for the interpretation and comparison 
with analytic results. In the following we review our results 1,2 from this 
perspective. 

QCD at finite quark chemical potential fi can be simulated with ordi- 
nary methods when /i is purely imaginary. If one were able to determine 
thermodynamic observables with infinite accuracy, standard complex anal- 
ysis arguments would guarantee that the result will be valid within the 
entire analytic domain, i.e. everywhere away from phase transitions. In 
practical numerical work one has to take into account two sources of er- 
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rors: first, the analytic form of the fitting function is not known a priori; 
second, even if it were so, one has to deal with numerical errors on the co- 
efficients. Cross checks among different analysis and guidance from models 
are thus most useful. Results from an imaginary \i have been obtained for 
the critical line of the two, three and two plus one flavor model 5 , as well 
as for four flavor . Thermodynamics results - order parameter, pressure, 
number density - were obtained for the four flavor model 2 . The pressure 
was defined by introducing an integral method at fixed temperature, and 
the mass dependence was estimated via the Maxwell relations n < 2 > 7 . 



2. The Critical Line, Chiral Symmetry and Confinement 

We used the exact results available for the critical line of the Gross Neveu 
model, and for Random Matrix Models in the appropriate universality class 
to show that a second order polynomial in T and fi approximates well the 
critical line over a large /i interval 1 . The nature of the chiral transition itself 
can be studied with a great accuracy : the correlation between < tyty > 
and Polyakov loop was observed at several values /Uj, either by studying 
the results as a function of /3, and by following the Monte Carlo histories 
directly at c (jJ>i) L ' 2 - 



3. The Hadronic Phase 

The grand canonical partition function of the Hadron Resonance Gas 
model 12 has a simple hyperbolic cosine behaviour. This can be framed 
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Figure 1. Hadronic Phase: (a) One Fourier coefficient fit to the particle number, show- 
ing that the Hadron Resonance Model is adequate to describe this data, (b) Compilation 
of the results for the chiral condensate and the particle number as a function of real 
chemical potential: the lines are cut in correspondence with fi c , showing the first order 
character of the phase transition (inferred from the chiral condensate) and the critical 
density. See our extended writeups for error analysis and details. 
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in our discussion of the phase diagram in the temperature-imaginary chem- 
ical potential plane which suggests to use Fourier analysis in this region, as 
observables are periodic and continuous there 1 . 

For observables which are even (O e ) or odd (0 o ) under fi — > — fi the 
analytic continuation to real chemical potential of the Fourier series read 
O e [o]{pi,N t ) — J2 n cosh[sinh](niV tJ ?V c ^/). In our Fourier analysis of 
the chiral condensate 1 and of the number density 2 - even and odd observ- 
ables, respectively - we limited ourselves to n — 0, 1, 2 and we assessed the 
validity of the fits via both the value of the x 2 /d.o.f. and the stability of 
a^p and a^) given by one and two cosine [sine] fits: we found that one 
cosine [sine] fit describes reasonably well the data up to T ~ 0.985T C (see 
Fig. la); further terms in the expansion did not modify much the value of 
the first coefficients and does not particularly improve the x 2 /d.o.f.. This 
means that our data are well approximated by the hadron resonance gas 
prediction AP oc (cosh(/is/T) — 1) in the broken phase up to T ~ 0.985T C . 
The analysis of the corrections requires better precision. 

The analytic continuation (Fig. lb) of any observable O is valid within 
the analyticity domain, i.e. till \i < ^ C (T), where /x c (T) has to be measured 
independently. The value of the analytic continuation of O at /i c , 0(/x c ), 
defines its critical value. When O is an order parameter which is zero in the 
quark gluon plasma phase, the calculation of 0(p, c ) allows the identification 
of the order of the phase transition: first, when 0(/j, c ) ^ 0, second, when 

o(/0 = o !' 2 . 

4. The Hot Phase 

The behaviour of the number density (Fig. 2a) approaches the lat- 
tice Stcphan-Boltzmann prediction, with some residual deviation. We 
parametrise the deviation from a free field behavior as 13,14 

AP(T, M ) = /(T, M )P/ ree (T, M ) (1) 

where Pf ree {T,fi) is the lattice free result for the pressure. For instance, in 
the discussion of Rcf. 14 

/(2» = 2(1 - 2a s /7r) (2) 

and the crucial point was that a s is /x dependent. 

We can search for such a non trivial prefactor f(T, /x) by taking the 
ratio between the numerical data and the lattice free field result n^ ree (fii) 
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at imaginary chemical potential: 

R{T^) = ^I\ (3) 

A non-trivial (i.e. not a constant) R(T,/j,j) would indicate a non-trivial 
/(?»• 

In Fig. 2b we plot R(T,fii) versus Hi/T: the results for 
T > 1.5T C seem consistent with a free lattice gas, with an fixed 
effective number of flavors N e / f {T)/-i = R{T): N e / f = 0.92 x 
4 for T = 3.5T C , and N e / f = 0.89 x 4 for T = 1.5T C . 
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5. Lattice vs. Models: 
Open Issues 

We give here a partial list of re- 
sults together with related ques- 
tions which calls for a more de- 
tailed understanding of the in- 
terrelations of analytic and nu- 
merical results. 

The critical line is similar to 
the one of model theories with 
the same global symmetries: 
T/T c 2 = 1 - 0.0021(2)(^/T) 2 
for four flavor, and similarly for 
two and three flavors 3 ' 5 ' 6 . Can 
we understand the flavor depen- 
dence and the (small) coefficient 
of (/z/T) 2 from these model? 

We have observed Tf lira \[i) = 

rpdeconfining^ Can e fJ ective 

Lagrangians 16 account for this 
observation ? 

In the hadronic phase 
AP ~ fc(l - cosh^s/T)). The 
corrections are not completely 
negligible, and can be used to 

estimate the magnitude of further terms in the expansion. Can these cor- 
rections be modeled in some simple, intuitive way? 



Effective Nf/4 T = 
Effective Nf/4 T = 



1.5 Tc > 
3.5 Tc ■ 



Figure 2. Hot Phase: (a) Particle number 
aproaching the free field behaviour (b) Ratio of 
the lattice results to the lattice free field: the 
deviation from free field can be described by an 
effective number of flavor N < 4 for T > 1 . 5T C 
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For T > 1.5T C the results are compatible with lattice Stefan Boltzmann 
with an active fixed number of flavor 4 x 0.92 for T=3.5 T c and 4 x 0.89 
for T — 2.5T C . These finite density corrections appear to leave unchanged 
the free field structure, even at moderately low temperatures. How can we 
understand this? 

Between T c and 1.5 T Cl can we interpret the deviations from this simple 
behaviour by use of a rigorous perturbative analysis 15 , and / or in the 
framework of a strongly interacting quark gluon plasma 17 ? 
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